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We demonstrate the existence of chaos in realistic models of two-field inflation. The 
chaotic motion takes place after the end of inflation, when the fields are free to oscillate 
and their motion is only lightly damped by the expansion of the universe. We then 
investigate whether the presence of chaos affects the predictions of two-field models, 
and show that chaos enhances the production of topological defects and renders the 
growth rate of the universe sensitively dependent upon the "initial" conditions at the 
beginning of the oscillatory era. 
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1 Introduction 



Inflationary cosmology and chaos theory are sites of intense activity within astrophysics 
and mathematics. As the equations governing the evolution of an inflationary universe 
are nonlinear, it is natural to ask whether the insight gained from studies of chaos in 
systems of differential equations sheds any light upon the behavior of the inflationary 
universe. In this paper we present prima facie evidence that the dynamics of realistic 
inflationary cosmological models based on two coupled scalar fields can be chaotic^] 

Chaotic dynamics in models consisting of a single (non-minimally coupled) scalar 
field in a Robertson-Walker background have been studied by several authors [^, ||, ||, |f|. 
Moreover, Cornish and Shellard || show that the familiar inflationary model of a 
minimally coupled field with a quadratic potential exhibits chaotic dynamics when 
spacetime has positive spatial curvature. 

In chaotic models with just one scalar field the scale factor must always be an 
intrinsic part of the chaotic system; thus, the timescale over which chaos becomes 
manifest is roughly defined by the lifetime of the universe. Cornish and Levin |7| point 
out that in two-field models the oscillations of the fields can be the primary source 
of chaos, so chaotic effects can become significant on comparatively short timescales. 
However, their analysis hinges on the dependence of the final state of the universe on 
the initial conditions, so their work does not fully realize the possibility that the time- 
scale associated with chaos in two-field inflation can be much less than the lifetime of 
the universe. Hence, while chaos has been found in several homogeneous inflationary 
models, it would typically only be detectable by an "observer" able to study either 
a large ensemble of universes, or the same universe during more than one cycle of 
expansion and contraction. A different situation arises in analyses of reheating, when 
the inflaton decays into a second scalar field ||, ^, 10 1. At least one of the fields is 
expanded in terms of its spatially dependent modes, so the resulting dynamical system 
is more complex than in the homogeneous case. However, small changes in the model 
parameters can cause massive variations in the occupation numbers associated with the 
modes, which is the qualitative behavior associated with chaos. 

The chaotic motion we examine is derived from the homogeneous equations of mo- 
tion and takes place in a well-defined era after the end of inflation. Consequently, 
the chaotic period we identify can, in principle, have effects which are observable from 
within the present-day universe. Furthermore, previously studied examples of chaos in 

inflationary models with chaotic dynamics, which are sensitively dependent upon their initial con- 
ditions, are not necessarily synonymous with chaotic inflationary models Jl|, which are recommended, 
in part, by their relative insensitivity to initial conditions. 
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homogeneous inflationary models typically depend upon non-zero spatial curvature, a 
cosmological constant, or a non-minimal coupling between the fields and the gravita- 
tional sector. All these terms are absent from our model, and the inflationary system 
we analyze has already received considerable attention because of its promising cos- 
mological predictions. The only additional assumption we need is that the parameter 
values and couplings between the scalar fields and the rest of the matter sector (which 
we do not include explicitly) permit a sufficiently long and lightly damped oscillatory 
period for the chaotic effects to become manifest. 

We begin by analyzing the evolution of the fields in the limiting case when the 
friction term arising from the expansion of the universe is absent from the equations 
of motion for the fields. We demonstrate that this system is chaotic, and that there 
is a minimum energy below which chaos does not occur. When the friction term is 
restored, we show that the energy changes slowly on a time scale defined by the typical 
oscillation period in the post-inflationary era. The chaotic properties of the system 
persist, but are now transient and cease to be apparent once the energy drops below 
the critical value needed for chaos. 

If inflation is the source of the fluctuations that seed structure formation, the "initial 
conditions" of the many post-inflationary horizon volumes accessible to a present-day 
observer cannot have been strictly identical. However, analyses of inflationary models 
usually make the tacit assumption that small differences in field values and velocities 
between post-inflation horizon volumes can be ignored when computing "zero order" 
quantities such as the scale factor. The validity of this assumption cannot be taken for 
granted in the presence of chaos. In particular, we consider whether chaotic motion 
enhances the production of topological defects, and if there is a well-defined relationship 
between the expansion of the universe and the decrease in density during the chaotic 
era. 



2 The Model and the Equations of Motion 

The equations of motion for inflationary models minimally coupled to Einstein gravity 
are well known, and for two scalar fields with a combined potential V((p, ip) they take 
the form 
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where a is the scale factor of the Robertson- Walker spacetime metric, m pl is the Planck 
mass, and dots denote differentiation with respect to time. The specific model we 
consider is defined by the potential 

m *) = (m 2 - ^V) 2 + + ^ v, (5) 

where we assume that 7 > 0. The inflationary properties of this potential have been 



widely studied [11]- [21]. Inflation is driven by the (p field: when <p is large, ip acquires a 
large positive effective mass, and the fields roll down the ip ~ groove in the potential. 
There are two distinct inflationary modes, which Copeland et al. [17] dub the inflaton 



dominated and vacuum dominated regimes. In the former case, the total energy density 
is much greater than M 4 , while during vacuum dominated expansion the energy density 
is of the order of M 4 . The degenerate minima located at (p = and ifi = ±M\/2A _1 / 4 
(separated by a barrier of height M 4 ) can lead to the formation of topological defects. 



The inflationary dynamics associated with this potential are described in 17], and 
sun 
where 



we summarize them here. The effective mass of the ip field becomes zero when 2 = <fif nst , 
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When 4> 2 first drops below (pf nst , ip is balanced on an unstable local maximum in the 
potential and becomes free to oscillate. Conversely, during inflation ip ~ and the 
potential reduces to the single field case, 

y(0)«M 4 + ^V. (7) 
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ceases when 



If inflation ends with <pr 3> (f>f nst , the M term can be ignored during inflation, which 



4. « (8) 
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and the slow roll conditions are no longer satisfied. If inflation continues until y - . - in _, . 
it ends quickly via the instability in the ip field. A realistic model in which inflation 
ends with the universe in the inflaton dominated stage must satisfy 

4^F— « 3x irr 3 , (9) 
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7 M 2 

7t > 8 *K- (10) 

5-5 x IP' 6 

m = -j= — m pl . (11) 

V87T 

The first inequality ensures that the M 4 contribution can be ignored, the second guar- 
antees that slow rolling is violated and inflation ends when 2 > cj) 2 nst , and the last gives 
COBE-normalized density perturbations .0 



3 The Frictionless Case 



We are concerned with the epoch in which the fields are weakly damped and free to 
oscillate. We begin our analysis by investigating the completely frictionless system. 
Physically, this is equivalent to assuming that the universe is not expanding and that 
the energy density is static. Mathematically, this means the equations of motion for 
the fields can be derived from a Hamiltonian. The equations of motion simplify further 
if we rescale the variables: 
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(12) 



The potential becomes 

V($,^) = M 4 

and the equations of motion are 
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where E is conserved. These equations are independent of M, so we fix M = 1 for 
convenience. If T = 0, we have a special case: the equations separate and our two 



dimensional system has two constants of the motion, ensuring that it is integrable [22|, 
and therefore not chaotic. To our knowledge, Hamiltonian chaos associated with the po- 
tential described by equation ( |i~3| ) has not been considered before. However, Matinyan 

2 The numerical value is taken from |r^| , and would be modified slightly by using the four year 
COBE result. However, the difference will not affect our work. 
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and Miiller |23|] discuss a system which lacks the <3? 2 /2 term in the potential but is 
otherwise equivalent to our model in the absence of friction. In general, Hamiltonian 
systems with terms like <1> 2 <I' 2 are non-integrable, and Steeb et al. summarize discus- 



sions of chaotic behavior when the potential consists solely of a <]? 2 ^ 2 term |24]. 
3.1 Applying the Toda-Brumer Test 



The Toda-Brumer test p5| , 26] checks for the presence of a growing, unstable solution, 
which is typically a prerequisite for chaos.0 The criterion for instability is the existence 
of and \E r for which 

\W"\\- d2Vd2V ( d2V ) 2 <0 (17) 
Inserting equation ( |l3| ) into equation ( |i~7|) gives 

I | = (r^ 2 + i)(3^ 2 -2) + r$ 2 (i-3r^ 2 ). (is) 



Setting equation (18) to zero and solving gives a parametric equation for the borders 



of the region within which ||V$^|| is negative, 

= gg-qpg+i) 

r(3r* 2 - 1) v ' 

Inserting this value of $ into V yields U(^f), the value of the potential on the 

boundary of the instability region, 

PW-fl-gUi-' 8 *'- 2 '' 1 *"'''. (20) 
K J \ 2 J 2T 3r^ 2 - 1 v ' 

Minimizing U(^>) subject to the constraint that $ is positive returns the minimum 
energy (as a function of T) needed for instability. When ty 2 = 2/3, $ = and U = 4/9, 
putting an upper bound on the minimum energy, which must be non-zero because 
IIKmJI * s positive at the minima of the potential. More precisely, solving U'^) = 
shows that U(^f) has extremal values at 

9 9r 9r v ; 

where we have ignored two roots for which is imaginary. A little analysis shows that 

when r < 9/2, the minimum energy for instability is 4/9, corresponding to <3? 2 = 0. For 

3 The Toda-Brumer criterion is equivalent to the Gaussian curvature of the potential being negative. 
However, even if ||V^|| is never negative chaos can still occur, so the Toda-Brumer test does not 
provide a sufficient condition for stability. 
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r > 9/2 the minimum energy for instability is found when \l/ 2 is given by equation (^l|), 
with the positive sign on the square root. As T — > oo, the minimal value approaches 
11/27, reproducing the result of Matinyan and Miiller in the absence of the <3? 2 /2 term. 

3.2 Lyapunov Exponents 

The Toda-Brumer test predicts instability above a critical energy, and other studies 
show that the coupling between the two fields usually leads to chaos. However, we wish 
to confirm this explicitly for the potential we consider, and to ascertain the parameter 
values which are likely to produce cosmologically significant chaos. 

We do this by considering the Lyapunov exponents |^7| . These are directly related 
to the fundamental definition of a chaotic system: sensitive dependence upon initial 
conditions. A system with an n dimensional phase space has n Lyapunov exponents. 
These measure the logarithm of the expansion of a small n- volume of phase space when 
it is propagated forwards in time by the equations of motion. In a chaotic system, 
this region will be stretched exponentially in at least one direction, and the signature 
of this is a positive maximal Lyapunov exponent. In Hamiltonian systems, the total 
phase space volume is conserved, leading to a spectrum of Lyapunov exponents that 
has the form Ai, A2, • • • , — A2, — Ai, [p^l - A corollary is that if a Hamiltonian system is 
not chaotic all of its Lyapunov exponents are zero. Since the Lyapunov exponents come 
in pairs with the same magnitude and opposite signs we can check that the values of 
the exponents we compute numerically are self-consistent. 

We calculate the Lyapunov exponents using the algorithm developed by Wolf et al. 
] which produces the full spectrum of n exponents at the expense of solving n 2 + n 
differential equations. In our case, n = 4, so the cost of calculating all the exponents 
is not prohibitive. However, Lyapunov exponents are defined by a limit as t —* 00 but 
the numerical integrations only run for a finite time. Hence our numerical estimates of 
the Lyapunov exponents are never precisely zero in non-chaotic cases. In practice, we 
continue the numerical integrations long enough to ensure that if Ai > 0.005 and Ai 
and IA4I differ by less than 1% the system is unambiguously chaotic. 

The maximal Lyapunov exponent as a function of E and T is displayed in fig. (|l]), 
while fig. (^) shows the region of parameter space for which chaos is observed. If 
r < 0.05, $ and are weakly coupled and the system is typically not chaotic. However, 
when E > 2 and T > 0.05 then we usually (but not always) find chaos. We very rarely 
found chaos with E < 1 and never observed chaos with an energy below the minimum 
value for instability obtained from the Toda-Brumer test. Furthermore, the boundary 
of the region in parameter space which gave rise to chaos is not smooth, so a small 
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change in parameter values could bring about a large change in the qualitative nature 
of the dynamics. For fixed E and Y only a subset of all the possible orbits in phase 
space will exhibit sensitive dependence, so the precise forms of figs (Q) and (|2|) depend 
on the chosen initial data. 



4 Chaotic Motion and Inflation 

Having demonstrated the existence of chaos in the Hamiltonian system derived from 
the two-field potential, equation (^), we now turn to the corresponding inflationary 
model based on the potential, equation (||). Because of the terms proportional to H in 
equations (||) and (|j), the energy density is strictly decreasing. If the energy changes 
slowly compared to the typical oscillation time of the fields, the friction term is a small 
perturbation to the Hamiltonian system. Since there is a lower bound on the energy 
at which the potential, equation (13), can exhibit chaotic motion, any chaos seen in 



the inflationary dynamics is transient, switching off once the energy drops below the 
critical value. 

In general, the oscillatory phase begins with ip ~ and <ft 2 ~ <fi 2 nst , although inflation 
may cease long before this condition is fulfilled. We define E inat , the (potential) energy 
at the beginning of the oscillatory epoch, 

y/Xm 2 \ „ r4 / H . 1 



E inst = V(<p inst , 0) = M \1 + = M ^1 + — J . (22) 

where F is given by equation (|l2|). The rate of change in the energy density is 

f^(f+?™H* #+n (23) 

where the second equality uses the equations of motion. During the oscillatory phase 
(cf> 2 + tp 2 ) ~ E. If we assume that M ~ m, the characteristic oscillation period of the 
fields when E < £7 inst and E ~ M 4 will be on the order of 1/M. Thus AE, the energy 
lost during a typical oscillation, is roughly 

^ ~ v^^. (24) 
E m p i 

Even if E inst is only a few times larger than M , and M ~ m ~ 10 _6 m pl , the fields will 
make tens of thousands of oscillations as the energy decreases from E inst to M 4 , when 
symmetry breaking occurs. This estimate is confirmed by direct numerical computa- 
tion. Thus the friction term is much smaller than the other terms in the equations of 
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motion, and the energy loss per cycle can easily be as low as 1 part in 10 5 . It is likely 
that the chaos seen in the Hamiltonian dynamics will survive the addition of such a 
small perturbation. 

We now turn our attention to the phenomenological aspects of chaos in the post- 
inflationary universe. We consider two ways in which the chaotic motion we have 
observed in the Hamiltonian dynamics could influence the observable properties of the 
universe: by enhancing the production of topological defects in the early universe, or 
by introducing stochastic variations into the relationship between length scales at the 
present epoch and the horizon size during inflation. 

The smallness of the perturbation caused by the friction term creates a practical 
difficulty Tracing the evolution of the fields through tens of thousands of oscillations 
is computationally expensive and, because the underlying system is chaotic, would 
require an unrealistic level of numerical precision. This problem is accentuated when 
we consider the evolution for many different choices of initial conditions. Thus we wish 
to make the friction more efficient while still ensuring that it is sub-dominant, so that 
the numerical calculations are tractable. We could do this by adding a coupling between 
the fields and radiation, in a manner analogous to that discussed by Albrecht et al. [29|. 
Alternatively, we can consider oscillatory motion at an energy scale considerably higher 
than that consistent with COBE so that H , which is proportional to the the square root 
of the energy density, is larger and the damping more efficient [ pCf . Both approaches 
reduce the length of the chaotic era, and therefore reduce the impact of the chaotic 
dynamics. We adopt the latter technique, as it does not alter the structure of the 
equations, and work with M = m = 1, A = 1, 7 = .5 and m pl = 5 x 10 3 m. We begin 
our integrations with E ~ E inst = 3M 4 . The typical duration of the oscillatory era is a 
few hundred oscillation times. The growth of the universe during the oscillatory era is 
not significantly changed by the scaling. 



4.1 Defect Formation 

Traditionally, only inflationary models which ended with a first order phase transition 
were thought to produce topological defects at the end of inflation. One of the new 
features of two-field inflation was that it can easily produce topological defects after 
a second order transition, provided the potential has more than one distinct vacuum, 
like the case we consider here. 

In fig. (||) we show the relationship between the sign of ip field at late times and the 
"initial" conditions that apply when E ~ E inst . We have to chosen to vary the initial 
velocities, while keeping the other parameters fixed. With 7 = 0.5 we see that even very 
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small changes in the initial conditions change the minimum of the potential picked out 
by ifi after symmetry breaking. Conversely, with 7 = large changes in the initial data 
are required to alter the final state, and there are sharp boundaries between regions 
of initial conditions space leading to different minima. This qualitative difference is 
attributable to the fact that with 7 = 0.5 we have transient chaos, whereas in the 
absence of friction the underlying dynamics are integrable with 7 = 0. 

The difference in initial energy between the different configurations with 7 = 0.5 
is 5E/E < 10 -4 , slightly more than that suggested by COBE, so the range of initial 
densities is not unreasonably large. The spread in energies for the initial conditions 
examined with 7 = is much wider, giving 5E/E < 5 x 10 -2 . 

We have chosen to work with parameters that lead to a much shorter oscillatory 
period than would be found with realistic values, so we underestimate the impact of 
the chaotic period. Furthermore, if <f> had rolled freely from the end of inflation its 
velocity would be much greater than the range of values we considered. This causes us 
to underestimate the initial energy and further reduces the duration of the chaotic era. 
On the other hand, we have ignored couplings between the scalar fields and the rest of 
the matter sector and spatial gradients in the fields. Incorporating these effects would 
provide additional friction which damps the oscillations of the <f> and ifj and therefore 
tends to suppress their chaotic motion. 

Because we do not include gradient terms, the plots in fig. (||) represent the symme- 
try breaking pattern for an ensemble of homogeneous universes (or horizon volumes), 
rather than for a single universe in which the field values vary significantly within a 
single horizon volume. Thus, we do not claim that the results obtained here are directly 
applicable to the rate of defect formation in the real universe. However, it is clear that 
two correlated horizon volumes which have not been in causal contact since the end 
of inflation are much less likely to remain correlated after symmetry breaking if the 
dynamics of the oscillatory phase permit transient chaos. Thus the chaotic properties 
of two-field inflation can be expected to enhance the level of defect production. 

These numerical results apply to a specific set of parameters, but in general the 
post-inflationary universe will be chaotic if the energy density at the beginning of the 
oscillatory period is high enough to produce chaos in the corresponding frictionless 
system. Thus transient chaos is a possibility if the parameters combine to give T > 
0.05 and the energy density is greater than than approximately 2M 4 at the end of 
inflation. If inflation ends with the universe in the inflaton dominated regime, cp 2 will 
be similar in size to E inst , equation (^), at the beginning of the oscillatory era and 
the energy requirement will almost certainly be satisfied. However, if inflation lasts 
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into the vacuum dominated regime the energy density at the end of inflation could be 
close to M 4 , which would not be high enough for chaos even if T, the rescaled coupling 
between the fields, was large. 



4.2 The Growth of the Scale Factor 

To make predictions from inflationary models we need to match comoving distances in 
the present universe to the comoving horizon size during inflation. The precise form of 
this relationship depends on the thermal history and effective equation of state that ap- 
plies after the end of inflation [31, 32]. Derivations of this relationship typically assume 



that the many post-inflation horizon volumes contained within the visible universe have 
grown by the same amount since inflation stopped. This assumption seems reasonable, 
as microwave background data suggests that the primordial universe is homogeneous 
to within 1 part in 10 -5 . However, a period of transient chaos can massively amplify 
small differences in initial conditions, so we now explicitly examine this assumption by 
studying the growth of the scale factor during the interval in which the energy density 
evolves from E inst to M 4 . We use the same sets of initial conditions considered in the 
previous subsection. 

In fig. (HI) we plot the magnitude of the scale factor, a, when E = M 4 , where a = 1 
when <p = 4> inst . In the chaotic case, an initial variance in the energy of less than .01% 
leads to differences of over 10% in the size of the universe at symmetry breaking. In the 
non-chaotic case the growth of the universe varies much less dramatically as the initial 
velocities of the fields are altered, and what difference we do see is attributable to the 
variation in initial energy corresponding to the different initial velocities. Moreover, in 
the non-chaotic case the growth of the universe between the start of the integration 
and symmetry breaking is a smooth function of the initial data, while for the chaotic 
case it is not. 

The chaotic case produces large variations in the scale factor because the precise 
mixture of kinetic and potential energy varies significantly with small changes in the 
initial conditions, as the fields settle into different meta-stable oscillatory modes. In 
turn, this alters the effective equation of state and the amount of growth that occurs 
as the energy density drops from E inst to M 4 . This is illustrated by fig. (||), where the 
evolution of tp is plotted for two slightly different sets of initial conditions. Obviously, 
the solutions diverge after a few oscillations, which is a consequence of the chaotic 
dynamics. Less dramatically, both solutions settle into a quasi-stable oscillatory mode, 
but the amplitude of the oscillations differs significantly between the two solutions. It is 
the latter difference that changes the average kinetic and potential energies. This effect 
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does not rely on the potential having more than one vacuum, which is a prerequisite 
for the formation of topological defects, but could occur in any two-field model with 
transient chaos. 

In practice, the variation in the scale factor due to chaotic motion occurring after 
inflation could introduce a stochastic term into estimates of the growth of the universe 
since the end of inflation. In turn, this would affect estimates of the power associated 
with the primordial fluctuations on different length scales, assuming that the fluctua- 
tions generated during inflation differ from the perfectly scale-free Harrison-Zel'dovich 
spectrum. The typical length over which the scale factor varies will be on the order of 
the horizon size at the end of inflation, which is vastly smaller than astrophysical length 
scales in the present universe. However, models which allow a short burst of inflation 
after the main period of inflation is over, such as thermal inflation |33| or more com- 



plicated two-field models [34|, could conceivably magnify the post-inflationary horizon 
volumes to the point where the stochastic properties of the scale factor would not be 
smeared out on astrophysical lengthscales. 



5 Discussion 

We have shown that the dynamics of the post-inflationary era in two-field inflation can 
be chaotic. We initially analyzed the evolution of the fields in the frictionless limit cor- 
responding to a static spacetime background. This reduced the equations of motion for 
the fields to a two-dimensional Hamiltonian system. By employing the Toda-Brumer 
test and calculating the Lyapunov exponents we demonstrated that the frictionless 
system is chaotic. This approach is free from many of the difficulties associated with 
chaos in general relativity, especially the co-ordinate dependence of the Lyapunov ex- 
ponents, as discussed by Rugh |3^]. We then demonstrated that the damping term 
induced by the expansion of the universe can be viewed as a small perturbation to 
the post-inflationary dynamics, and that the most important qualitative property of 
chaos - sensitive dependence on initial conditions - persists when the fields evolve in 
an expanding universe. The chaotic era is transient and lasts from the beginning of 
the oscillatory period, which commences when the effective mass squared of the ip field 
becomes negative, until the energy density, which is strictly decreasing, drops below 
the critical value needed to sustain chaos. 

Our analysis is based on a specific potential, which we chose because its inflationary 
properties had already received widespread attention, and because it allowed us to 
investigate the formation of topological defects after the end of inflation. From the 
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perspective of the idealized, frictionless dynamics chaos is endemic in two-field models 
with couplings like ^(j) 2 Tp 2 . Thus we expect that the type of chaos described here could 
arise in any inflationary model with two coupled fields which undergo lightly damped 
oscillations. 

In models of reheating and particle production based on parametric resonance the 
inflaton field is coupled to another boson, usually a scalar, and the resulting system 
is similar (and sometimes identical) to the two-field model considered in this paper. 
Refs ||, |l0| show that particle production during reheating can depend delicately 
upon the model parameters, which is the qualitative behavior associated with chaos. 
These calculations expand the field(s) in terms of their spatially dependent modes, so 
the underlying dynamical system is more complex than the homogeneous case discussed 
here. However, the link between the chaos we have demonstrated in the homogeneous 
equations of motion and the phenomenology of reheating is extremely intriguing. 

In practice there are several ways in which chaos might be suppressed in a realistic 
model of two-field inflation; for instance, the parameter values could ensure that the 
oscillatory period does not begin until the energy density is too low to permit chaos, 
or the couplings between the scalar fields and other parts of the matter sector may 
ensure that the oscillations are always strongly damped. Conversely, our numerical 
integrations were performed with an artificially large damping term, and the parameter 
values we chose ensured that the oscillatory period did not begin until the energy density 
was close to the height of the barrier between the two vacua in the potential. Restoring 
a realistic damping term or increasing the energy density at which oscillations begin 
would make the chaotic period last much longer than it does in the case considered 
here. 

We have not included spatial variations in the fields, or couplings between the 
scalar fields and other forms of matter. Thus we cannot determine whether a chaotic 
period necessarily leads to detectable effects, and if these effects are useful to the model 
builder, or in conflict with observational evidence. However, our work shows that it 
is feasible for chaotic effects to influence the observable properties of two-field models. 
Thus, approximations which do not preserve the chaotic properties of the full system 
could lead to spurious conclusions about the phenomenology of two-field inflation. It 
will obviously be important to investigate these issues in more detail. 

Inflationary models based on a single scalar field in a spatially flat background do 
not possess enough degrees of freedom to become chaotic. Since spatial curvature is 
negligible after inflation, a chaotic era at the end of inflation requires two or more scalar 
fields. Thus any distinctive observable effects associated with the period of transient 
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chaos we have identified could provide a mechanism for testing a broad class of two- 
field models, in the same way that constraints on the scalar and tensor perturbation 
spectra test a large group of models with a single slowly rolling field |36|]. Furthermore, 
chaotic effects are often associated with fractal structures. Due to the complex nonlin- 
ear dynamics associated with defect formation, cosmic string networks have a fractal 
geometry |}7], [38] , while stochastic inflationary inflationary models have a fractal struc- 
ture at scales much larger than a typical horizon []39| . Thus we speculate that a chaotic 
era at the end of inflation may give rise to fractal structure, in which case two-field 
models of inflation could make testable predictions about the fractal properties of the 
universe at large scales. 

We believe that the scenario analyzed here is the most interesting example of chaotic 
dynamics in an inflationary model to date. We have found chaos in a model that 
has already been widely studied because of its promising phenomenology, and we did 
not need to include non-minimal couplings or non-zero spatial curvature to ensure 
the presence of chaos. Because of the sensitive dependence of chaotic systems, small 
differences in the initial conditions which apply at the beginning of the oscillatory 
epoch can have a significant impact on the subsequent evolution, affecting "zero order" 
quantities such as the scale factor and the minimum of the potential that the fields 
evolve towards. Moreover, because the chaotic era occurs after the end of inflation its 
effects might, in principle, be detectable by human observers. 
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Figure 1: Maximal Lyapunov exponent, Ai as a function of energy and T. The specific 
initial data chosen for the integrations was 

$ = 1, ^ = i/y/2 and $ = 0.64$. The 
magnitude of the velocities is then fixed by E, and their signs were chosen to be positive. 



16 



0.5 1 1.5 2 2.5 3 3.5 

r 



Figure 2: Region of parameter space for which chaos was detected (white points), with 
the same parameters as used in fig. ([!]), for a grid of 840 x 840 equally spaced points. 
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0(0) 

Figure 3: The minima of the potential reached at late times, as a function of initial 
conditions (for 420 x 420 equally spaced points), with 0(0) = inst and "0(0) = 0. White 
points evolve towards the vacuum with -0 > 0. The upper plot shows the results for 
7 = 0.5, and those for 7 = are shown in the lower plot. Note the much larger range 
of initial data considered in the 7 = case. We have chosen 0(0) negative so that is 
initially rolling "downhill". 



18 




Figure 4: The scale factor at the moment the energy density becomes equal to M 4 
and symmetry breaking occurs is plotted as a function of initial conditions, for 7 = 0.5 
(left) and 7 = (right). The scale factor, a, is normalized to be unity at the beginning 
of the integration. The initial data are the same as those used for fig. (||). 
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Figure 5: The evolution of ij) is plotted for two different initial conditions. In both cases 
the parameter values are the same as those used for the chaotic solutions in fig. (||), 
with the specific initial initial conditions V = -02 and <j> = —0.01(1 + 1/419) (upper) 
and 4> = —0.01(1 + 2/419) (lower). The particular numerical values have been chosen 
to reproduce specific points in fig. (|j). 
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